The shape of a small universe: signatures in the cosmic microwave background 
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We consider the most general parametrization of flat topologically compact universes, complement- 
ing the work of Scannapieco, Levin and Silk to include non-trivial shapes. We find that modifications 
in shape of the fundamental domain will lead to distinct signatures in the anisotropy of the cosmic 
microwave radiation. We make a preliminary assessment of the effect on three statistics: the angu- 
lar power spectrum, the distribution of identified "circles" on the surface of last scattering and the 
correlation function of antipodal points. 



PACS Numbers : 98.80. Cq, 98.70.Vc, 98.80.Hw 

With the dramatic improvement of cosmological obser- 
vations in the past few years, efforts are underway to 
precisely characterize the nature and structure of space- 
time. One property which has generated much interest 
is the possibility that we live on a manifold which has a 
compact, non-trivial topology |Q. A multitude of tests 
and observables have been proposed which would allow 
us to pin down the topological class of the universe we 
live in. Under the assumption that we live in a smooth 
Friedman- Robertson- Walker universe with perturbations 
set up adiabatically in the early universe, measurements 
of the Cosmic Microwave Background anisotropy seem to 
indicate that our universe is flat || . This greatly reduces 
the set of models one needs to look at and simplifies the 
boundary conditions of the wave functions one needs to 
consider. Scannapieco, Levin and Silk || have identified 
complete sets of functions for each topology for the case 
where the identification vectors are orthogonal. 

In the context of Kaluza-Klein theories, it has recently 
been pointed out in |^||| that non-orthogonality of the 
identification vectors will lead to interesting effects in the 
mass differences of Kaluza-Klcin states. It was shown in 
H that the mass gap of the graviton will become depen- 
dent on what he calls the shape parameter, i.e. the angle 
between identification vectors. While he considered the 
case of a compact 2-dimensional space rather than the 
standard 3-dimensional space we live in, we shall see what 
happens if we apply such a transformation to a compact- 
ified version of our 3-dimcnsional space in the case of no 
extra dimensions. 

Without loss of generality let us consider the simplest, 
flat, non trivial topology- the hypertorus in the classifi- 
cation of || . This will allow us to outline the formal- 
ism which can be used in the general case. Consider a 
3- manifold with coordinates ^=(^1,^2,^3) on which any 
function, /(£i,£2,£3) satisfies the following property: 



/(&+27r,&,&) = /(£i,6 + 27r,&) 
= + 2tt) = /(£i, £2, &) 



(1) 



We can think of this manifold as tiled by cubes with 
sides of length 2tt. Alternatively we can think of it as 
simply consisting of a cube with opposite faces identified; 
for convenience let us call it a fundamental domain. We 
can now rescale each direction in such a way as to have 



different identification scales in each direction by defining 
a matrix L: 



L = diag( J Ri,i? 2 ,i?3) 



(2) 



and a new set of coordinates: \ — L • £. It will be useful 
to divide this transformation into two: one that is simply 
the product of the identity by some overall scale factor, R 
and another transformation which is diagonal with unit 
determinant. The first transformation corresponds to an 
overall volume transformation, modifying the volume of 
the fundamental domain by a factor of R 3 . The sec- 
ond transformation induces anisotropic scale differences 
in the sides of the fundamental domain (changing it into a 
paralliliped). A coordinate invariant definition orders the 
eigenvalues of this second transformation, Ai > A2 > A3 
and defines o~\ — A2/A1 and o~2 = A3/A1. 

We now wish to consider an additional linear transform 
that changes the directions of the identification vectors. 
Naturally we are not interested in overall rotations of 
all vectors, given that ultimately any observable we con- 
struct should be rotationally invariant. Consider a set 
of 3 vectors: si, s*2, S3. Given that we wish to factor 
out overall rotations, we can align s\ with £1 , have s*2 lie 
on the (^1,^2) plane and S3 will have an arbitrary direc- 
tion. We can define a shape transformation matrix, S, 
such that its three columns are the three vectors, s,-. A 
convenient parametrization is the following: 



1 cos OL\ cos «2 
■j: I sinai sin CV2 cos 0:3 
sin a.2 sin 0:3 



(3) 



Note 

that the the prefactor, uj — (sinai sincv2 sin as) -1 / 3 , en- 
sures that det S = 1 . This parametrization is such that 
s 1 • S2 = to 2 cos a.\ , s 1 ■ S3 = uj 2 cos ct2 and S2 • S3 = 
w 2 (cosai cos 02 + sinai sin 0:2 cos 03). The geometrical 
interpretation of this transformation is simple: we are 
changing the angles of the corner of our fundamental do- 
main while preserving its volume. We can now consider 
a new set of coordinates x — SL • £. The identification 
vectors are not orthogonal anymore. 

Note that we are allowed to perform these transfor- 
mations without jeopardizing the suitability of the man- 
ifold as a cosmological model. We start off with a locally 
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homogeneous and isotropic space time with non-trivial 
identifications which are permitted. We then simply ap- 
ply a linear transformation which preserves the homo- 
geneity and isotropy of the manifold as well as differ- 
entiability. Indeed if we consider a generalization to a 
flat manifold in N dimensions with analogous topological 
constraints, we are allowed to perform any non-singular 
linear transformation and the resulting manifold will still 
have the same differentiability class. Such a transforma- 
tion will be characterized in terms of N rescaling param- 
eters and N(N — l)/2 shape parameters. Furthermore, 
these transformations clearly cannot change the topology 
class of the manifold. In the example we are considering, 
the topology remains that of a hypertorus and the geom- 
etry remains the same: zero curvature. 

The plane wave solutions to the Laplaces equations will 
be modified. A complete set of orthogonal solutions is 
given by ipn(0 — exp(— m-£) with n=2ir(ni, 712,713). This 
can be mapped directly onto a complete set of orthogonal 
solutions on the transformed space: 4>%{x) = exp(— ik ■ 
•?) = where the last equality implies k = L _1 S _1 n. 

As above the wavenumbers are quantized, but the effect 
of S is to mix up the quantization constraints through 
the shape angles. The explicit form of the wavenumbcr 
for this topology is 

wR\k\ cota 2 

= TL\ — 77,2 COt OLl + 71 3 (COt «1 COt «3 
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In this note we wish to focus on the effect of shape. We 
will restrict ourselves to Ri = i? 2 = R3 = R and two 
simple, one paramater families of shape transformations: 

A: We fix a 2 = 03 = ir/2 and vary ol\ = a e [0, n/2]. 
This essentially corresponds to collapsing one of the 
identification vectors of the fundamental domain on 
the plane defined by the other two. 



B: We fix si • s 2 
a e [0,7r/2]. 



§2 ' S3 = S3 • Si = cos a and vary 



An hypothesis was put forward in Q that the mini- 
mum traversable distance of the space defines the scale at 
which observable effects appear in statistics of the CMB 
such as the angular power spectrum or the geometry and 
frequency of ghost images. Although our shape transfor- 
mation has unit determinant, the minimum traversable 
distance (d m i n ) may be greater or less than the one in 
the cubic domain. For the two families we above we find 
that 
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for cos a G [0, 0.5] 



dmin — Rlo 



\/2(l — cos a) for cos a € [0.5, 1] 



(5) 
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FIG. 1. The minimum traversible distance for the two fam- 
ilys of shape transformations we are considering, as a function 
of the cosine of the angle. The solid (dashed) line corresponds 
to type A (B). 

and in Fig |l| we plot this quantity for the two families 
we are considering. We find that, indeed, d m i n does vary 
as a function of shape parameter. Although there is a 
slight increase for cos a < 0.5, both parametrizations A 
and B lead to a dramatic decrease in d m in for cos a ~ 
1. The question we now wish to address is, how will 
this be refelected in some of the statistics that have been 
proposed to quantify topology from the CMB. 

Let us first focus on the angular power spectrum of 
the CMB. In a flat universe, the anisotropy in the CMB, 
AT (6) (where is the unit vector pointing in a given 
direction in the sky), is dominated by the Sachs- Wolfe 
term, i.e. 



AT(0) = -i$(^,0d* 



(6) 



where $(77, x) is the gravitational potential, 77* is the 
time of last scattering, and is the comoving distance 
to the surface of last scattering. We will be working 
in Fourier space where /(x) = /(fc) exp(— i/c • x). 
Assuming statistical homogeneity and isotropy we have 
that ($*(r) ,k)<f>*(r)0,k')) = V{k)5 %% ,. Defining a lm = 
J dQT(6)Y[ m (6) we have that 

(a* im a e < m >) = (4:TT) 2 i l '- l J2v(k)j e (kd*)je,(kd*) 
k 

xY lm {k)Yl, m ,{k) (7) 

In general, for arbitrary topology the covariance matrix 
of the ai m will not be diagonal. The fact that there are 
preferred orientations in space (defined by the identifica- 
tion vectors) will induce correlations in the a£ m between 
different I and m which depend on the orientation of the 
Ytmifi) basis. In this letter, as in we wish to look at 
the overall impact on a rotationally invariant measure of 
the angular power spectrum. We shall consider 
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£ \a £m \ 2 = (4n)J2nk)jUkd*) (8) 
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FIG. 2. The angular power spectrum for a selection of three 
values of the shape parameter for the two families of trans- 
formations, A and B. The angle a takes values 7r/2 (dotted), 
7r/4 (long-dashed) and 7r/8 (solid). 

For the infinite volume case this converges to the stan- 
dard angular power spectrum. 

In Figure @we plot angular power spectra for a — 7r/2, 
7r/4, and 7T/8. We are considering a universe with an 
identification scale R — d*. For the cubic domain the 
features identified by || are manifest, i.e. one finds a 
supression of power for low I and a set of oscillatory fea- 
tures out to large I. Changing the value of a will shift the 
oscillatory features around but can have a curious effect 
at the lowest multipoles: the solid lines (corresponding to 
a = 7r/8) show an increase in power due to the fact that 
the shape transformation can lead to the regeneration of 
small wavenumber modes in the integral. Indeed we find 
this to be a general trend. 

All real space statistics rely, to some extent, on the rep- 
etition of patterns on the surface of last scattering. The 
surface of last scattering will intersect points in space 
which are identified. This leads to perfect correlations 
between certain points (or pixels) of a map of the CMB. 
The more points which are identified, the more signifi- 
cant a detection of topology will be. We can attempt 
to quantify the effect of shape by estimating how many 
points (or structure, such as for example circles) will be 
identified for a given value of a in either of the transfor- 
mations, A and B. Two points x and x' lie on the surface 
of last scattering if they satisfy \x\ 2 — \x'\ 2 — d 2 and can 
be identified with one another if there exist a trio of in- 
tegers, (m, 712,^.3) such that x! — x + n\Si + 712S2 + n 2 S2- 
Combining these two sets of equations one finds that a 
given trio of integers satisfies this condition if 

(n\ +n\+n\ + 2nxn 2 si ■ s 2 + 2n 2 n 3 s 2 ■ s 3 (9) 
+2n 3 m§s-h) 1/2 /[h-(hxss)] 1/3 <^. 



For example, for an orthogonal trio of Si and R = 2d* 
one has 3 pairs of identified structures. This condition 
for transformation A is: 

{n\ + n\ + n\ + 2n Y n 2 cos a)/ (sin a) 1/3 < ^ (10) 

R 

and for transformation B is: 
[n\ + n\ + n 2 + 2cosa(nin 2 + n 2 n 3 + n^ni)} 1 ^ 2 / 

(l-cosa)y(l + 2cosa)] 1/3 < (11) 

R 

As cos a — * 1, more roots are possible, leading to stronger 
correlations between points on the surface of last scatter- 
ing. 

There is in fact a simple way of understanding how the 
number of roots (or matched structures) increases with 
cos a. Let us first consider transformation A. The direc- 
tion of the minimum traversable distance is given (with 
cos a > 0.5) by si — S2. As we increase cos a, we de- 
crease dmin so that the number of identified structures 
will be given by the integer part of 2d*/d m in- Given that 
linicoso,^! d min = 0, there will be more matched struc- 
tures with increasing cos a. For transformation B, the 
reasoning is the same except there are now three possi- 
ble directions: Si — s 2 ,S2 ~ S3 and S3 — s\. We then have 
that the number of matched structures is given by the 
integer part of 6d*/d m i n . 

Let us now focus on two real space statistics that have 
been advocated for identifying topology in the CMB: 
matching circles in the sky byjTJ] and the correlation func- 
tion of antipodal points of ||jT In [Q Cornish, Spergel, 
and Starkman have proposed matching the circles which 
arise from the intersection of the surface of last scat- 
tering with the identified boundaries of the fundamental 
domain. One way of thinking about this is to consider a 
given identification vector s. One can then construct two 
perpendicular planes at each end and find the lines of 
their intersection with the sphere of radius d* and center 
at the midpoint of the vector. For example, for a cubic 
domain of size R < 2d* one expects at least 3 pairs of cir- 
cles with centers on the 3 axes of identification. Pairs of 
circles will intersect (with each other) if R < \/2d*, and 
each circle will subtend an angle 9q = ~ — arcsin -M- on 
the sky. For a cubic domain only pairs of big circles will 
intersect. Let us now consider a fundamental domain 
with a non-trivial shape parameter. The identification 
vectors will not be orthogonal and it will be possible to 
have pairs of small circles that intersect. The condition is 
that a < 9o ■ This is then a signature which could distin- 
guish between cubic and non-cubic domains. More gen- 
erally, if one identifies more than one pair of circles and 
the vectors joining their centers are not perpendicular, 
then there is evidence for non-trivial shape. Naturally 
one would have to distinguish it from the hexagonal fun- 
damental domains, where the identification vectors are 
at 60°. 
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FIG. 3. A one dimensional slice of the antipodal correlation 
function for a selection of three values of the shape param- 
eter for the two families of transformations, A and B. The 
angle a takes values n/2 (solid), n/4 (short-dash) and n/8 
(long-dash). 



Levin et al. in || have suggested that a good statistic 
for identifying non-trivial topology is the correlation of 
antipodal points on the surface of last scattering: 

C A (0) = (AT(-6)AT(§)) oc ^ V(k) cos(2<2*k • 0). (12) 



The patterns which appear will reflect the geometry of 
the fundamental domains. Again one might expect that 
a non-cubic domain may lead to a different signature; 
for example it has been shown in [0 that a hexagonal 
fundamental domain will have a very different signature 
than a cubic one. As an illustration of this statistic for 
non-trivial shapes we have calculated Ca{9) for the one 
dimensional slice 9 = (cos <fi, sin <fi, 0). We can cleary see 
the differences in Fig. g where we plot Ca{4>) f° r the two 
families of shape transformations. 

In this report we have taken a first look at how the 
shape of the fundamental domain of a flat, topologically 
non-trivial universe might affect a few different statis- 
tics of the CMB. There clearly seem to be observable 
signatures which will constrain the values of the shape 
parameters. We will now finish with a few comments. 
To begin with, it is clear that when it is claimed that one 
can use the cosmic microwave background to constrain 
the topology of the universe, it is in fact meant that one 
is constraining certain geometric properties which arise 
when the universe is topologically compact. For exam- 
ple, if the identification scales are infinite, then the shape 
transformations are irrelevant. When choosing the topol- 
ogy of the universe (from the six possibilities discussed 
in H ) , one must also include the 6 parameters which de- 
scribe the geometry and scale of the fundamental domain. 



One can then proceed to apply the variety of statistics 
that have been advocated. It is possible, however, to 
incorporate the parameters describing the fundamental 
domain into the standard likelihood methods which have 
been used to constrain cosmological parameters. The 
standard parameters pi for Gaussian models, should be 
extended to include another six parameters (R, o~\, a<x, 
ot\, a.2, 6*3) for topology. 

Topologically non-trivial universes are no longer statis- 
tically isotropic. One must choose an overall orientation 
of the fundamental domain. This complicates consider- 
ably the evaluation of the likelihood function for a given 
topology primarily because one cannot use the radical 
compression methods || . One must work with the full co- 
variance matrix in real space, C{9, 6), which is no longer 
simply a function of 9 ■ 9. This is a daunting computa- 
tional task: not only is one extending parameter space by 
another six factors, but one must also marginalize over 
all possible orientations in space. Furthermore, what was 
once a one-dimensional radial integral over wavenumber 
becomes a three dimensional integral which only picks 
up contributions from the selected modes on the Fourier 
grid. If one is to correctly constrain the properties of the 
fundamental domain of our universe, it is essential to de- 
velop novel techniques which can incorporate these vari- 
ous subtleties into likelihood methods. One should point 
out, however, the impressive increase in speed of the algo- 
rithms for performing the various steps in maximum like- 
lihood estimation. Hopefully this means that the chal- 
lenge of constraining topologically non-trivial universes 
will be met and that it will be possible to do so with the 
data which will come from the MAP satellite. 
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